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1. Introduction 

Coulomb gauge is a suitable choice for investigating the confinement phenomenon, as in this 
gauge the Gribov-Zwanziger scenario becomes distinctly prominent: the temporal component of 
the gluon propagator provides the long range confining potential, whereas the spatial propagator is 
infrared suppressed ||l|]. On the other hand, traditional studies of the heavy quark sector of quantum 
chromodynamics (QCD) mainly use phenomenologically motivated potentials in the place of the 
Yang-Mills sector. In this context, it is appropriate to investigate the relationship between the 
nonperturbative scale associated with confinement and the Yang-Mills sector of the theory. 

This talk reviews results obtained in the heavy quark limit of Coulomb gauge QCD [Q, ^. Af- 
ter making an expansion in the heavy quark mass and restricting to the leading order, we consider 
the heavy quark propagator and the homogeneous Bethe-Salpeter equation for quark-antiquark sys- 
tems. With the further truncation to exclude pure Yang -Mills vertices (but retaining nonperturbative 
dressed propagators), we show that the rainbow-ladder approximation is exact in this case and we 
establish a connection between the temporal gluon propagator and the external physical scale (the 
string tension), at least within the leading order truncation. In the second part of the talk, we inves- 
tigate the (full nonperturbative) four-point quark-antiquark Green's functions. We present exact, 
analytic solutions, and show that the physical poles of the Green's function expUcitly separate from 
the possible unphysical ones. Moreover, we find that the physical poles of the Bethe-Salpeter equa- 
tion are contained within the singularities of the Green's function. These results will hopefully be 
useful in the future investigations of phenomenological models for mesons and baryons (see, for 
example, Ref. [Q] for a numerical analysis of the inhomogeneous Bethe-Salpeter equation). 

2. Quark propagator in the heavy mass limit 

Let us start by considering the explicit quark contribution to the full QCD generating func- 
tional: 

= j ^<J>exp|; j d^xq„{x) ifOo + ifD-m ^^^/^wj 

X exp I J J d'^x [Xa {x)qa (x) + ^„ {x)Xa i^)] + i^YM^ ■ (2.1) 

In the above, i^<I> denotes the integration over all fields present, qa is the quark field, qa the con- 
jugate antiquark field, and Xa,Xa the corresponding sources. The common index a, j3 . . . denotes 
the color, spin and flavor indices. The Dirac / matrices satisfy {y^,y'^} = 2g^^, with the metric 
gf^^ = diag(l, — 1). The structure constants of the SU{Nc) group are denoted with f"^'^, and the 
Hermitian generators satisfy [T'',T^] = if'^^'T'^ and are normalized via Tr(r",r^) = 5"^/2. ^ym 
represents the Yang-Mills part of the action and 

D^) = ^Q-lgT''o"{x), D = V + igr'A"{x), (2.2) 

are the temporal and spatial components of the covariant derivative in the fundamental color repre- 
sentation (A and a refer to the spatial and temporal components of the gluon field, respectively). 
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We now decompose the full quark field according to the heavy quark transformation 

'l + f 



qa{x) = e-"'"" [h{x) +H{x)]^ , ha{x) = e" 



-q{x) 



, Ha{x)=e 

a 



imxQ 



(2.3) 



(similarly for the antiquark field), where the two components h and H are introduced with the help 
of the spinor projectors (1 it 'f')/2. This is a particular case of the heavy quark transform, adopted 
from the Heavy Quark Effective Theory [HQET] [||], which turns out to be useful in Coulomb 
gauge. 



After inserting the quark fields, decomposed according to Eq. (2.3), into the generating func 



tional Eq. (2.1), we integrate out the //-fields and make an expansion in the inverse of the heavy 
quark mass (in the following, we will adopt the standard terminology and denote it simply "mass 
expansion"). At leading order, the generating functional reduces to: 



'[la] = y"^<I>exp|; j d\ha{x)[idox + gT"a"{x)]^php{^ 

xexpiij d'x [e-'"'">Xa{^)ha{x) + e"^''ha{x)Xa{x)] + i^YM \ + ^(1/m) , (2.4) 



where we have replaced the covariant derivative Dq with its explicit expression. In the above 
expression, we notice that as a result of the heavy quark transformation, at leading order in the 
mass expansion the quark interacts only with the temporal gluon, whereas the spatial component 
is suppressed. Also, note the absence of the Dirac 7 structure, as a result of the multiplication 
with the projectors (1 ± 7°)/2 (physically, this implies that the spin degree of freedom decouples 
from the system). A further important point is that in Eq. ( |2.4| ) we have kept the full quark and 
antiquark sources, as opposed to HQET, where the sources corresponding to the large components 
h are used. This means that at leading order in the mass expansion we are allowed to use the full 
apparatus of the functional formalism, and hence derive the full Dyson-Schwinger equations in 
Coulomb gauge QCD, while replacing the corresponding propagators and vertices by their leading 
order expressions. 

In Coulomb gauge QCD (without the mass expansion), the quark gap equation for the proper 
two-point function is given by (see Ref. for a complete derivation and notation): 

=rfl,{k) +ld(o[r^^^^pik,-(o,(o-k)w^^p,{(o)^^,^^^^ 

where dco = d(o/{2nY. In order to solve this equation, we first need to specify the tree-level 
quark proper two-point function and the components of the tree-level quark-gluon vertex. They are 



derived directly from the generating functional Eq. ( |2.4| ) and are given by: 

4?«/3 (^) = ['^o - m] + ^ ( 1/m) , (2.6) 
r£i«/3(^i,^2,^3) = [gr^p + ^(1/m) . (2.7) 
Note that the spatial component of the quark-gluon vertex is of order ^(1/m), as shall be explained 



shortly below. Further, the nonperturbative temporal gluon propagator entering Eq. (2.5) has the 
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form [0]: 



(2.8) 



Following lattice results, which signal that that the gluon dressing function is largely independent 
of energy, we assume that Dca is a function of the three-momentum. Also, lattice investigations 
indicate that Daa is infrared divergent and behaves like 1/P for vanishing k^ [^] (however, we will 
need the explicit form of this function only at the end of the calculation). 

Finally, the last input is provided by the Slavnov-Taylor identity, which furnishes a relation 
between the two- and three-point functions of the theory. This is derived from the invariance of 
the QCD action under a time-dependent Gauss-BRST transform and in Coulomb gauge reads 
{ki+k2 + k^ = 0): 



+ 



aS 



Sp 

^qqSp{-h) 



(2.9) 



In the above, qo is an arbitrary energy injection scale (arising from the noncovariance of Coulomb 
gauge [0]), Tec is the ghost proper two-point function, and Tq-ccq and Tq-ccq are ghost-quark kernels 
associated with the time-dependent Gauss-BRST transform. Since in the generating functional, 
Eq. (2^), the tree-level spatial quark-gluon vertex appears at ^(1/m), by making the further trun- 



cation to neglect the pure Yang-Mills vertices, it follows from the Dyson-Schwinger equation for 
the spatial quark-gluon vertex that the fully dressed F^^a ~ ff{l/m) (and hence it is neglected). 
Further, the ghost-gluon vertices involve pure Yang-Mills vertices and hence are also truncated out. 
Thus, in our truncation scheme and at leading order in the mass expansion, the Slavnov-Taylor 
identity takes the simple form: 



'^3^qqaalii''uh,k3) ='i^qqas{h) 



sp 



aS 



^qsp{-h) + ^{l/m) 



(2.10) 



as a solution of Eq. (|2.5|), combined with Eq. ( |2. 10| ) 



Collecting the above results, we find the following expression for the heavy quark propagator. 



qqap 



(k) 



[ko — m — J^r + 1£ 

with the constant (implicitly regularized, as indicated by the index "r"): 



^r = \g^C,j- 



«2 



+ ^(l/m), 



(2.11) 



(2.12) 



where rfw = d^a>/(27c)^ and Cf = {N^ — l)/2Nc- The gap equation, Eq. ( p3| ), has been solved 
under the assumption that the temporal integral has been performed first, with the spatial integral 



regularized and finite. The solution, Eq. (2.11), is then inserted into the Slavnov-Taylor identity, 
Eq. ( Z.lCp , and we find that the temporal quark-gluon vertex remains nonperturbatively bare: 



^qaapikuk2M) = [gT"]„p + ^{y/m). 



(2.13) 
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Note that the quark propagator, Eq. (2.11), possesses a single pole in the complex plane, 
in contrast with the standard QCD propagator, which has a double covariant pole. Hence, we 
need to derive the quark and antiquark propagators separately, with the corresponding Feynman 
prescriptions. Examining the closed quark loops (virtual quark-antiquark pairs connected by a 
primitive vertex), we find that they vanish due to the energy integration over two quark propagators 
with the same Feynman prescription, and this implies that the theory is quenched in the heavy mass 
limit: 



[ko-m — + 1£\ [k{)+ pQ — m — + 1£\ 



0. (2.14) 



Further, note that the propagator Eq. ( |2. 11| ) is diagonal in the outer product of the fundamental 
color, flavor and spinor spaces, and this exhibits the decoupling of the spin from the heavy quark 
system. Finally, it is important to emphasize that the position of the pole has no physical meaning, 
since not the propagator itself, but the bound state of a quark and an antiquark is physical. The fact 
that the poles in the quark propagator are shifted to infinity once the regularization is removed sim- 
ply means that an infinite energy is needed to create a single quark from the vacuum. If a hadronic 
state is considered, only the relative energy (derived from the homogeneous Bethe-Salpeter equa- 



tion) is required to describe the system, and in this case the singularities appearing in Eq. ( |2.11 ) 



cancel. Similar types of cancellation appear in the solutions of the four-point quark-antiquark 
Green's function (see also the discussion from section^. 
For the antiquark propagator we obtain: 



[ko + m — J^r + l^] 



Wg^ap{k)= ,, Z . .^, +^(l/m), (2.15) 



and the corresponding temporal antiquark-gluon vertex is given by: 

r^5cTa/3(^i>^2,^3) = -[gT"]p„ + ^{l/m). (2.16) 

In the above, notice the Feynman prescription of the propagator, as well as the sign of the loop 
correction, J^,-. As shall be discussed in the next section, this apparently minor modification will 
play an important role in the interpretation of the solutions of the Bethe-Salpeter equation for 
quark-antiquark states as bound state/confining solutions. 



3. Homogeneous Bethe-Salpeter equation 

Let us now consider the full homogeneous Bethe-Salpeter equation for quark-antiquark bound 
states [Q] (see also Fig. [T]): 

nP-^P)aP = dkK^p.sMk-P) [Wgg{k+)r{k;P)Wgg{k^)]^s. (3.1) 

In the above, the momenta of the quarks are given by p+ = p + ^P, = p — {I — (similarly 
for k±), and is the momentum sharing fraction (note that the solutions here and in the next section 
turn out to be independent of i^, just as in the covariant case [p^). P represents the 4-momentum 
of the bound state (assuming that a solution exists), F is the Bethe-Salpeter vertex function for the 
particular bound state and K is the Bethe-Salpeter kernel, which still needs to be specified. 
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Figure 1: Homogeneous Bethe-Salpeter equation for quark-antiquark bound states. Internal propagators are 
fully dressed and solid lines represent the quark propagator The box represents the Bethe-Salpeter kernel K 
and filled blobs represent the Bethe-Salpeter vertex function Y with the (external) bound state leg given by a 
dashed line. 



It is well-known that the Bethe-Salpeter kernel K and the quark self-energy are related via the 



axialvector Ward-Takahashi identity [|11[]. Under truncation, it was seen in the last section that the 
heavy quark self-energy reduces to the rainbow truncated form; the corresponding truncation for 
the kernel is the ladder approximation. In this has been explicitly derived. The Bethe-Salpeter 
kernel in the heavy quark-antiquark system under the truncation considered here is thus: 



We now insert the nonperturbative results for the propagators and vertices, Eqs. ( |2.11| , [2. 13 



n|,|[T6|) and the expression Eq. ( |2.8[ ) for the temporal gluon propagator. After explicitly identi- 



fying the antiquark contribution, i.e. Wqq{k_) = —Wqq{—k-), we perform the temporal integration 
over the quark and antiquark propagators, which now leads to (unlike Eq. ( 2.14j )): 



^"^^ (3.3) 



Further, we insert the expression Eq. (2.12) for J^^ and after Fourier transforming to coordinate 



space we find the following solution for the bound state energy of the quark-antiquark system: 



UP- 



Cp — e' '^Cm 



+ ^(l/m). (3.4) 



In the above. Cm is an (unknown) color factor assigned to the Bethe-Salpeter vertex F, which has 
to be yet identified: 

[T"Y{x)r\^=CMTafi{x). (3.5) 

Because the total color charge of the system is conserved and vanishing [p^], a quark cannot 
exist as an asymptotic state. Hence, the bound state energy Pq of a quark-antiquark system can 
be either infinite, such that the system is not allowed, or linearly rising, if the system is confined. 
If the temporal gluon propagator is more infrared divergent than l/|w|, we find that in order to 
ensure the convergence of the spatial integral. Cm must be equal to Cp- This immediately leads to 
the condition 

Tay{x) = 5ayT{x), (3.6) 

which implies that the Bethe-Salpeter equation can only have a finite solution for color singlet 
states. Further, if we assume that in the infrared Daa =X/gP (as indicated by the lattice [j8|]). 
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where X is a combination of constants, then from Eq. (3.4) we find 



Po = o\x\ 



871 



\x\ + &{l/m) 



(3.7) 



This result shows that there is a direct connection between the physical string tension a and the 
nonperturbative Yang-Mills sector of QCD, at least under the truncation scheme employed here, 
which corresponds to a color singlet bound state of a quark and an antiquark, and otherwise the 
system has infinite energy. 



4. Four-point quark-antiquark Green's functions 

Similar to the gap equation, we derive the the full Dyson-Schwinger equation for the proper 
(IPI) four-point quark-antiquark Green's function in configuration space (see Ref. [||] for notation 
and technical details of the functional derivation): 



<i'laxi<lrziqT:wi1nt>= [gfT"]ap J dy5{x-y) 



><iq^iqy-jO;^><ip;^ipS> 



< 



lXpxlXS><iqsiqvtl(^e><lp£lpK> 



<iq^^.iqvioj[><iZyiXe><i(]eiQr]t^(^S><^P^^Py> 
<iq^„iqv iO^^><iJyiXx> <iqx i^yz^^^i > < «Pp ipy > 



<iXpxiXK><iqKi^rzi(^x > <^Px 'Pv > <qTwiqmiOyia^ > <ipf, ip" > 



+ 
+ 



<apx^XS><m^(lr]tl<ye ><WeWK> 



<iq^„iqy-JO"^lO'i><iplip;> 



<iXpxiXK><iqK"lrziQtwiQnt(y'l><ip^^ipy> 



< 



iXpx^XKXiqK^lrz^qxi^ntXiqTwiflvif^iiXiXviXxXip^ip 



> 



<iXpxiXK><iqKi(i7ziqTwm><iXx^Xs><i(is^<iritiol><iplipy> 



< 



iXpxiXv><iqviqiiiq^„iqr]t><ix^iXK><iqKiqyzi(^x><^Px^P} 



> 



+ 



+ . 



< 



iXpx^XK><iqK^qjziOx><iPxipy 



> 



<iq^„iq^iay><iX^iXS><tqsiqritl<ye><iPsiPy> 



<apx^xs><miqr]ti(ye><WeWv> <iqTwiqni(^v><ix^iXic><iqKiqrzi(y;^><ip^ipy>^ j- 

(4.1) 



where the dots represent the A vertex terms (which will be truncated out in our scheme), and we 



have already replaced the tree-level temporal quark-gluon vertex with its expression Eq. ( |2. 13| ). 
The above equation is diagrammatically represented in Fig. ^ 

We now proceed by applying our truncation scheme at leading order in the mass expansion. 
Since we shall consider the flavor non-singlet Green's function in the s-channel (the quark and the 
antiquark are regarded as two distinct flavors, but with equal masses), the diagrams (a), (c) and (i) 
of Fig. are excluded. In the diagram (b) (crossed ladder type exchange diagram), we insert as 



before the appropriate propagators Eqs. (2.11, 2.15) and vertices Eqs. (2.13, 2.16). The resulting 
energy integral is similar to the integral (2.14) and vanishes, just like the higher order contributions 
to the kernel of the homogeneous Bethe-Salpeter equation. Turning to the diagram (d), we see 
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(f) (g) (h) (i) (j) 



Figure 2: Diagrammatic representation of the Dyson-Schwinger equation for the IPI 4-point quark- 
antiquark Green's function. Blobs represent dressed proper (IPI) 4-point vertex, solid lines represent the 
quark propagator, springs denote either spatial (A) or temporal (a) gluon propagator and cross denotes the 
tree level quark-gluon vertex. Internal propagators and IPI vertices are fully dressed. 



that this contribution involves a quark-two gluon vertex. From the corresponding Slavnov-Taylor 
identity (derived explicitly in Ref. |^]), it can be seen that this vertex is zero and hence the diagram 
(d) vanishes. 

At this stage, we adopt the following strategy: discard for the moment the diagrams (f) and (g), 
which include the IPI four-point quark-antiquark Green's function, and the diagram (e), containing 
a four-quark-gluon vertex, solve the equation with the remaining terms (the diagram (h) and the 
rainbow-ladder term (j)), and with the obtained solution return to the diagrams (f), (g) and (e), and 



show that they cancel (and hence our assumption is consistent). In this case, Eq. ( |4. 1| ) reduces to 
the Dyson-Schwinger equation for the IPI four-point quark Green's function in the 5-channel, in 
the ladder-approximation (shown diagrammatically in Fig. 

^alz7i{PhP2,P3,P4) = 

-J ^fqaap(Pl^-Pl-(^^^)^q^PsiPl+0))^^aSiiiPl+(^^P4,-Pl-P4-(0] 



Wf^i-(0)W^i{py+P4 + (0 



- J d (Or^^^^^p{pi,-pi - (0,(0)Wgqpy{pi + (0)r'ClrT^ipi + (0,P2- (0,P3,P4)W^^^^^ 

xr|,^^y(« -P2,P2, -co)Ca(-«)- (4-2) 



As before, we identify the antiquark component of Eq. (4^) (lower line of Fig. |3|) and insert 



the expressions Eqs. (2.11, 2.15), for the quark and antiquark propagators, along with the vertices 
Eqs. (2.13, 2.16| ) and the definition Eq. ( ^ ) for the temporal gluon propagator. We further make 
the assumption that r''^\pi,p2,P3,P4) = r(^^(Po;pi + P4), where Po = p^ + P2^ which allows us 



to separate the three-momentum and energy integrals. The energy integrals are similar to Eq. (3.3) 
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- — CO 




ra+p,+p. 




X 



Figure 3: Truncated Dyson-Schwinger equation for the IPI 4-point Green's function in the i-channel. 
Same conventions as in Fig. 2 apply. 



and can be carried out. Making the following color decomposition for the function F^^^ : 



^(4) 
ayTt] 



(4) 



(4) 
2 ' 



(4.3) 



where Pl^^ and T 



^(4) 



2 are scalar functions and Fourier transforming back to coordinate space, we find 
the following solution for the IPI quark-antiquark Green's function: 

2 

I 



2Nr 



o2 



X <^ 5ay5. 



Po - 2J^r + l^Waaix) + lie 
(PO - 2^r)Nc{Nl - 2) + Ig^NcCpWaaix) 



Po - IJr - Ig^CpWaaix) + lie 

where x = |x| is the separation associated with the momentum p\ + p\. 



xy 



(4.4) 



As promised, with the solution Eq. ( |4.4| ) for the IPI Green's function, we now return to the 
diagrams (f), (g) and (e) and show that they do not contribute to the final result. To see this, we first 
consider the diagram (g) and notice that the energy dependence of the internal four-point function 
can be written as 



0. 



(4.6) 



where X is a combination of constants, ?i = 1,2 and m = 0, 1. Then the energy integral takes the 
form 

2+n Y 

Clearly, this integral is a generalization of Eq. ( |2.14[ ) and this vanishes, just as for the loop correc- 
tions in the kernel of the Bethe-Salpeter equation from the previous section and the diagram (b) 
from above. An identical calculation for the diagram (f), recalling that the lower line corresponds 
to an antiquark propagator, leads us to the fact that this integral is also zero. Finally, turning to 
the diagram (e), containing the four quark-gluon vertex, we notice that the perturbative series of 
this diagram coincides with the ladder resummation of the diagrams (f) and (g), which we have 
found to be vanishing, and hence this diagram is also zero (even though the five-point interaction 
vertex itself does not vanish - see also Ref. [||] for a detailed diagrammatic analysis). In turn, this 



implies that our original assumption is correct and the solution Eq. (4.4) is valid at every order in 
perturbation theory. Marginally, we note that the fact that the five-point function is finite relates 
to the existence of three-quark bound states in the Faddeev equation, in the ladder approximation 
Pll. 
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Figure 4: Relation between the IPI (dark blob) and amputated (shaded blob) 4-point Greens function for 
the quark-antiquark system. Internal propagators and IPI vertices are fully dressed. 




Figure 5: Truncated Dyson-Schwinger equation for the fully amputated quark-antiquark 4-point Green's 
function in the i-channel. 



Let us now consider the Dyson-Schwinger equation for the fully amputated four-point quark- 
antiquark Green's function in the s-channel, which we denote G'^^\ This is related to the IPI 
function r'"*) via the Legendre transform (see Fig. ^). This study is motivated by the fact that this 
equation reduces (under truncation) to the inhomogeneous ladder Bethe-Salpeter equation, from 
which the homogeneous Bethe-Salpeter equation of the previous section is derived. In the follow- 
ing, we will derive the solutions of this equation, analyze the positions of the poles, and explic- 
itly verify that the physical solutions coincide with the bound state solutions of the homogeneous 
Bethe-Salpeter equation. 



Starting with Eq. (4.2) (in coordinate space) for the proper function, we replace the IPI func- 



tion r^^^ with the amputated function G*^^' according to Fig. ^ and cut the external quark legs. The 
resulting equation for G*^^' reads (see also Fig. ^: 



G%,r,{p+,P--M,k-) = w^*'„{p-k) [r| 



i arj 



qqa 



zy 



dq 



MO) 



Wqq{ci^ 



Wqqiq-)r 



qqa 



{p - q)GKli;.n {q+ ,q-;k+,k^). (4.7) 



In the above, for the quark momenta we use the same conventions as in the Bethe-Salpeter equation, 
Eq. (3.1). Again, we replace the heavy quark and antiquark propagators and vertices with the 



expressions Eqs. ( |2. 1 1| , [2.15| , [2.13| , [2. 16| ), perform the energy integration, and Fourier transform 
back to coordinate space. Similar to the proper four-point function, we make a color decomposition 
of the function G'^^^ : 



^(4) 



(4) 



(4) 



(4.8) 



{G\ and Gj are scalar functions), and after using the Fierz identity to sort out the color factors. 



^(4) 
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we obtain the final result for the function G^*^ : 



G 



(4) 



iP0-2.yr)Waa{x) 



2 Po - + ijN^^oo (x) + lie 



' Pq - IJy - ig^Cp^aa {x) + lie 
is the total energy of the qq state. 



Nr 



(4.9) 



where Pq = 7^+ 

Analyzing the structure of the four-point functions, we first notice that even though the results, 
Eqs. ( |4.4| , [4.9| ), are derived under truncation, the denominator structure of the IPI and amputated 
Green's functions are identical in both color channels, and furthermore, the physical and nonphys- 
ical singularities have disentangled automatically. Using the form Eq. ( ^!^ for the temporal gluon 
propagator, the denominator factor of the color singlet channel in either Eq. ( pf ) or Eq. ( p3| ) can 
be rewritten in the form 

W rfo3D (03) r 
Jr or- L J 

In this expression we recognize the bound state (infrared confining) energy Pq ,-es {x) = o \x\ (fur- 
ther assuming that D(j(j(w^) ~ 1/w^), similar to Eq. (3.4) for the homogeneous Bethe-Salpeter 



equation in the color-singlet channel. Hence, we have found an explicit analytical dependence 
of the four-point Green's function on the qq bound state energy resulting from the homogeneous 
Bethe-Salpeter equation. 

Turning to the overall denominator factors of Eqs. ( p!4l , i.e. the denominator factor not 
specific to the color-singlet channel, we again insert the explicit form of the temporal gluon propa- 
gator and arrive at the following result: 



da)D(ja{(o) 

^2 



Cf + 



1 



(4.11) 



This factor does not appear in the homogeneous Bethe-Salpeter equation; it is part of the normaliza- 
tion and, similar to the quark propagator, represents an unphysical pole which is shifted to infinity 
when the (implicit) regularization is removed. 



5. Conclusions 

In this talk, we have discussed the Dyson-Schwinger and Bethe-Salpeter equations for quark- 
antiquark systems in Coulomb gauge, at leading order in the heavy quark mass expansion and 
with the truncation to include only the (nonperturbative) temporal gluon propagator. Under this 
truncation, the rainbow approximation to the quark gap equation is exact, as is the corresponding 
ladder approximation to the homogeneous Bethe-Salpeter equation. The only physical solution 
corresponds to confinement, i.e. only color singlet meson states have finite energy (and hence are 
physically allowed), and otherwise the system has infinite energy. Incidentally, these results are 



supported by recent Dyson-Schwinger studies in Coulomb gauge at leading order []14|]. 

Turning to the four-point quark-antiquark Green's functions, we have presented analytic solu- 
tions for both the proper and amputated Green's functions. The two functions have the same de- 
nominator structures, and the physical and nonphysical singularities disentangle, the physical poles 
coinciding with the bound state solutions obtained for the homogeneous Bethe-Salpeter equation. 



11 



Nonperturbative study of the four-point heavy quark Green's functions 



Carina Popovici 



Acknowledgments 

CP. has been supported by the Deutscher Akademischer Austausch Dienst (DAAD), the Helm- 
holtz International Center for FAIR within the LOEWE program of the State of Hesse, and the 
Helmholtz Young Investigator Group No. VH-NG-332. P.W. and H.R. have been supported by the 
Deutsche Forschungsgemeinschaft (DFG) under contracts no. DFG-Re856/6-2,3. C.R thanks the 
organizers, in particular C. Aguilar and D. Binosi, for the support. 

References 

[1] V. N. Gribov, Quantization of non-Abelian gauge theories, Nucl. Phys. B 139 (1978) 1; D. Zwanziger, 
Lattice Coulomb Hamiltonian and static color-Coulomb field, Nucl. Phys. B 485, 185 (1997); 
D. Zwanziger, Renormalization in the Coulomb gauge and order parameter for confinement in QCD, 
Nucl. Phys. B 518 (1998) 237. 

[2] C. Popovici, P. Watson and H. Reinhardt, Coulomb gauge confinement in the heavy quark limit, Phys. 
Rev. D 81 (2010) 105011. 

[3] C. Popovici, P. Watson and H. Reinhardt, Higher order heavy quark Green' s functions in Coulomb 
gauge, Phys. Rev D 83 (2011) 125018 [arXiv: 1103.4786 [hep-ph]]. 

[4] M. Blank and A. Krassnigg, Matrix algorithms for solving (in)homogeneous bound state equations, 
Comput. Phys. Commun. 182 (2011) 1391 [arXiv: 1009. 1535 [hep-ph]]. 

[5] M. Neubert, Heavy quark symmetry, Phys. Rept. 245 (1994) 259. 

[6] C. Popovici, P. Watson and H. Reinhardt, Quarks in Coulomb gauge perturbation theory, Phys. Rev. 
D 79, 045006 (2009). 

[7] P. Watson and H. Reinhardt, Two-Point Functions of Coulomb Gauge Yang-Mills Theory, Phys. Rev. 
D 77, 025030 (2008) [arXiv:0709.3963 [hep-th]]. 

[8] M. Quandt, G. Burgio, S. Chimchinda and H. Reinhardt, Coulomb gauge ghost propagator and the 
Coulomb potential, PoS CONFENfEMENT8, 066 (2008). 

[9] P. Watson and H. Reinhardt, Slavnov-Taylor identities in Coulomb gauge Yang-Mills theory, Eur. 
Phys. J. C 65 (2010) 567. 

[10] R. Alkofer, P. Watson and H. Weigel, Mesons in a Poincare covariant Bethe-Salpeter approach, Phys. 
Rev D 65 (2002) 094026 [hep-ph/0202053]. 

[11] S. L. Adler and A. C. Davis, Chiral Symmetry Breaking in Coulomb Gauge QCD, Nucl. Phys. B 244 
(1984) 469. 

[12] H. Reinhardt, P. Watson, Resolving temporal Gribov copies in Coulomb gauge Yang-Mills theory, 
Phys. Rev. D79 (2009) 045013. 

[13] C. Popovici, P. Watson and H. Reinhardt, Three-quark confinement potential from the Faddeev 
equation, Phys. Rev. D 83 (2011) 025013 [arXiv: 1010.4254 [hep-ph]]. 

[14] P. Watson and H. Reinhardt, Leading order infrared quantum chromodynamics in Coulomb gauge, 
arXiv:l 11 1.6078 [hep-ph]. 



12 



